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ABSTRACT 

Gravitational waves (GWs) from supermassive binary black hole (BBH) inspiral s are potentially 
powerful standard sirens (the GW analog to standard candles) ljSchutznl986l 12002) . Because these 
systems are well-modeled, the space-based GW observatory LISA will be able to measure the luminos- 
ity distance (but not the redshift) to some distant massive BBH systems with 1-10% accuracy. This 
accuracy is largely limited by pointing error: GW sources generally are poorly localized on the sky. 
Localizing the binary independently (e.g., through association with an electromagnetic counterpart) 
greatly reduces this positional error. An electromagnetic counterpart may also allow determination of 
the event's redshift. In this case, BBH coalescence would constitute an extremely precise (better than 
1%) standard candle visible to high redshift. In practice, gravitational lensing degrades this precision, 
though the candle remains precise enough to provide useful information about the distance-redshift 
relation. Even if very rare, these GW standard sirens would complement, and increase confidence in, 
other standard candles. 

Subject headings: black hole physics — gravitation — gravitational waves — galaxies: nuclei — cosmology: 
observations — cosmology: theory — gravitational lensing 



1. INTRODUCTION 

One of the major challenges to cosmology for the 
foreseeable future is to understand "dark energy" , the 
mysterious component responsible for the apparent ac- 
celerating expansion of our Universe llRiess et alJll998t 
[ Perlmutter et all Il999t iTonrv et all l2003t iKnop et all 
2003). Dark energy can be parameterized by its con- 
tribution to the universe's energy density, Six, and 
its equation-of-state ratio, w{z). Of particular in- 
terest will be measurements that probe w(z), testing 
whether the dark energy is a true cosmological constant 
[w(z) = —1] or wh ether it arises, for example, from an 
evolving field (e.g. iCaldwell. Dave, fc Steinhardtl Il998t 



Type la SNe are excellent standard candles, with a 
(calibrated) peak brightness thought to be known to 
about 15%. A possible objection to SNe as standard 
candles is the absence of a solid theoretical underpin- 
ning. Of particular concern is the possibility of evolution 
in S N brightnesses, leading to unknow n systematic er- 
rors ijDrell. Loredo. fc Wassermanll2000|) . In this article, 
we discuss a completely independent standard candle: 
the gravitational-wave (GW) driven inspiral of massive 
binary black holes (BBHs). As GW detections can be 
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One of our best observational probes of the dark en- 
ergy is the distance-redshift relation, which maps the 
expansion history of the universe. Much of our knowl- 
edge of this relation comes from observations of dis- 
tant Type la supernovae (SNe). These SNe serve as 
standard candles: their observed intensity can be cali- 
brated to tell us their luminosi ty distance, Dj, JP hillit>s 
Cini IRiess. Press, fc KirshnerllS IWang et a flT20n3L 
As the redshift of a SN (or its host) can also be mea- 
sured, each SN puts a point on the distance-redshift 
curve. Future surveys (e.g., SuperN ova/ Acceleration 
Probe 1 , Large- aperture Synoptic Survey Telescope 2 ) are 
expected to measure thousands of Type la SNe, mapping 
the distance-redshift curve over a large span of redshift 
with good statistical significance. 



http: //snap . lbl . gov 
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thought of as aural rather than optical (Hughes 2003), 
a more appropriate term for a GW standard candle is a 
"standard siren". 3 Because BBH systems are relatively 
simple and well modeled (at least in the early "inspi- 
ral" phase of their coalescence), the GWs they generate 
determine the source's luminosity distance with high ac- 
curacy: typically 5Dl/Dl ~ 1-10%, with most of the 
uncertainty ar ising from correlations with pointing errors 
(Hughes 2002). BBH merger events will follow the merg- 
ers of galaxies and pregalactic struc tures at high redshift 
ijVolonteri. Haardt. fc Madaul2003|) . Though the merger 
rate is poorly understood, LISA is expected to mea- 
sure at least several events over its mission, especially 
as it is sensitive to these waves to enormous distances 
()Richstonelll99l lHaehneltj|l99^ . 

Since GWs do not provide the redshift of the source, 
BBH GW measurements alone do not probe the distance- 
redshift relation. However, as first noted by Bernard 
Schutz, should some kind of "electromagnetic" (EM) 
counterpart to a BBH GW event be identified, the situ- 

3 We thank Sterl Phinney and Sean Carroll for suggesting this 
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ation changes drastically l)Schutdll986l 120021. First, by 
determining the source position, many correlations which 
set the distance error are broken. The error then drops 
immensely — below 0.5-1% in many cases. Second, a 
counterpart could determine the source's redshift. A 
BBH GW source coupled with an EM counterpart could 
therefore constitute an exceedingly good standard siren 4 . 
We comment at this point that, to date, there has not 
been a great deal of careful analysis regarding the nature 
of EM counterparts which may accompany a GW event. 
We discuss briefly some ideas that have been presented 
to date regarding the form that counterparts may take in 
SjHland 351 We hope that the promise of this high quality 
candle will motivate additional thinking on this issue. 

In practice, gravitational lensing will limit the qual- 
ity of this candle. GWs are lensed by inter- 
vening matter exactly as electromagnetic waves are 
lensed (lMarkoviclll993t [Wang. Stebbins. & Turnerlll996t 



iTakahashi fc Nakamura 12003(1 . As the waves propagate 
through our inhomogeneous universe, they are magnified 
(or demagnified) , inducing some error in our inferred lu- 
minosity distance to the source. As we discuss in ^ the 
distribution of errors is such that a BBH candle will most 
likely be comparable in quality to a Type la SN standard 
candle. However — and we strongly emphasize this point 
— the BBH candle will have entirely different systemat- 
ica from SNe. Concordance between the two types of 
measurement could thus alleviate concerns about evolu- 
tionary effects in Type la SNe, and greatly increase one's 
confidence in all standard candles. 

2. DISTANCE DETERMINATION WITH BBH GWS AND 

LISA 

Massive BBH coalescences are among the most lumi- 
nous events in the universe. That luminosity (peaked at 
~ 10 57 erg/sec) is radiated in GWs, which couple very 
weakly to matter. The planned space-based GW detec- 
tor LISA (the Laser Interferometer Space Antenna) will 
be sensitive to these BBH waves in the frequency band 
(10" 5 -10- 4 ) Hz < / < 0.1 Hz, making possible measure- 
ments fr om binaries with t otal masses mi + 7712 ~ 10 3 - 
1O 6 M {D anzmanr] 11998(1 out t o redshifts of at least 
z ~ 5-10 and possibly beyond l|Hughesl l2002t IVecchiol 
2004). In this section we discuss how LISA measure- 
ments determine the distance to a source, summarizing 
our model of the waveform and LISA's sensitivity and 
response, and discussing the measurement precision we 
expect from measuring merging black hole populations. 

2.1. Merging black hole GWs 

For this paper, the most interesting epoch of BBH co- 
alescence is the inspiral, when the binary's members are 
widely separated and slowly spiral together due to back- 
reaction from GW emission. The GWs from this epoch 
are well modeled using the post-Newtonian approxima- 
tion, roughly speaking an expansion in inverse separation 

4 It may also be possible to use the distributions of observed 
binaries for cosmology, obviating the need for an EM cou nterpart 
IChernoff He FinrlHwa IFinrilT99a I Wang fc TurneillTgSll . Unless 
the event rate is much higher than currently expected, however, 
the statistical errors associated with these distributions suggest 
that these methods will not achieve accuracy sufficient to measure 
properties of the dark energy equation of state, our primary focus. 
Certainly other cosmologically interesting measurements could be 
made. 



of a binary's members; see iBlanchetl l)2002fl and refer- 
ences therein for more detailed discussion. We will not 
discuss waves from the merger (in which the holes come 
into contact, forming a single body), nor from the ring- 
down (the final, simple stage of the ringdown, in which 
the merged binary is well-modeled as a single, distorted 
black hole) , as they do not substantially impact distance 
determination. 

Inspiral GWs encode the luminosity distance to a bi- 
nary, its position on the sky, its orientation, and informa- 
tion about c ertain com binations of masses an d spins; see 
lArun et al\ (J2004) and IBlancheTeT al. ( 2004) for up-to- 
date discussion and details. The inspiral does not encode 
a source's cosmological redshift. Redshift is instead en- 
tangled with the binary's evolution. For example, the 
masses (mi,m2) impact orbital evolution as timescales 
(Gmx / c 3 , Gm 2 / c 3 ) . These timescales redshift, so the 
measured masses redshift: a binary with masses (mi , 7712) 
at redshift z is indistinguishible from a local binary with 
masses [(1 + z)m\, (1 + z)ma] (modulo amplitude). This 
reflects the fact that general relativity has no absolute 
scale. 

In a reference frame centered on the solar system's 
barycenter, the strongest harmonic of the inspiral GWs 
two polarizations has the form 



hy = 



2M 5 J 3 [7rf(t)]y 3 r 

D L 

AMl /3 [irf{t)] 2 / 3 (L-n) 



1 + (L ■ n f cos[$(i)] , (1) 
sin[$(t)] . (2) 



The mass parameter M z = (1 + z)(mim,2) 3 / 5 /( m i + 
ni2) 1 ' 5 is the binary's redshifted "chirp mass", so called 
because it largely sets the rate at which the binary's 
members spiral towards one another, determining the 
"chirp" of the orbital frequency. The phase $(t) de- 
pends on intrinsic binary parameters — the ma sses and 
spins of its members (e.g.. iPoisson fc Willlll995l) . It de- 
pends particularly strongly on M. z \ as a consequence, 
phase coherent measurements of the wa yes will determine 
the chirp mass with great pre cision l|Finn fc Chernofj 
119981 iCutler fc Flanagan! ^994) . The wave frequency 
fit) — (l/2-?v)d$>/dt. The unit vector n points from 
the center of the barycenter frame to the system, and 
hence defines its position on the sky; L points along the 
binary's orbital angular momentum, and hence defines 
its orientation. Notice that the luminosity distance Dl 
appears in combination with these two angular factors. 
Determining Dl thus requires fixing these angles. As we 
now discuss, LISA is able to do so by virtue of its orbital 
motion. 

2.2. Merger GWs as measured by LISA 

The LISA antenna consists of three spacecraft, ar- 
ranged in orbits about the Sun such that they form an 
equilateral triangle (roughly; the armlengths are in gen- 
eral not equal, and in fact oscillate — albeit with periods 
much longer than that of the GWs we aim to observe). 
This triangle "rolls" as the spacecraft move through their 
individual orbits, preserving the triangular formation. 
The centroid of the constellation shares Earth's orbit, 
lagging by 20°, so that it takes 1 year to orbit the Sun. 
Figure ^ shows a schematic of the orbital configuration. 
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Fig. 1. — Illustration of the LISA antenna's orbit. The con- 
stellation "rolls" as its centroid orbits the sun, completing one full 
revolution for each orbit. 



See|D anzmanr] l)1998(l for detailed discussion of the LISA 
mission and its orbital configuration. 

At least in the low frequency limit (/ < c/L, where 
L is armlength), LISA can very usefully be regarded as 
two GW detectors: the time varying armlength data, 
(SLi, 8L2, SL3) can be used to synthesize outputs for two 
equivalent "L" -shaped detectors, with 90° arms. These 
"equivalent dete ctors" are rota ted by 45° with respect to 
one another; see lCutlerl l)1998j) for details and derivation 
of this viewpoint. The datastream sm of the these two 
equivalent detectors is given by a weighted sum of the 
two GW polarizations, plus noise: 

*Ui(t) = \ PiTn(*)M*)+^i*nM*)J +ni,n(t) • (3) 

The prefactor y/3/2 in this expression enters when con- 
verting the "real" interferometer response to that of the 
synthesized equivalent detectors. The antenna functions 
F^~j* depend on the orientation and position of the 
source relative to the antenna. Because of the antenna's 
orbital motion, the position and orientation of the source 
relative to the antenna is continually changing. The mo- 
tion of the detector thus modulates the measured signal; 
the exact nature of the modulation depends upon the 
position and orientation of the source. We write the re- 
sponse functions as time dependent functions to reflect 
this modulation. Note also that the waveform phasing 
is modified in an important manner by the antenna's 
motion — the orbital motion causes freque ncy as well 
as amplitude modulation. See iCutlerl (|1008f) for further 
discussion. 

We take the noises in the equivalent detectors, niji(t), 
to be uncorrelated, Gaussian random processes, with the 
same RMS values: 

(mmi) = ; (n?) - (n 2 u ) . (4) 

In all of our analysis, we use t he sam e noise model as 
that used by iBarack fc Cut led <|2004D [their Eqs. (48)- 
(54)]. In our calculations, it is necessary to introduce a 
low frequency cutoff — a frequency at which the sensi- 
tivity to GWs rapidly degrades. This cutoff has impor- 
tant implications for determining which binaries LISA 
can measure: the frequency support of a binary's GW 
spectrum is inversely proportional to its mass. In other 
words, more massive binaries will radiate at lower fre- 
quencies than less massive binaries. The low frequency 
cutoff thus determines the maximum binary black hole 



mass accessible to LISA measurements. It also deter- 
mines the amount of time for which a binary's waves are 
in band: a binary that may only be in band for a few days 
when /i ow = 10 -4 Hz may be in band for many months 
when /i ow = 3x 1CT 5 Hz. 

Unless stated otherwise, we have set /i ow = 1CU 4 Hz for 
the results we present here. This is a somewhat conser- 
vative choice; some members of the LISA mission design 
community (particularly P. Bender) argue that LISA 
should have good sensitivity down to frequencies / ~ 
1CU 5 Hz. Accordingly, we have put /i ow = 3 x 10~ 5 Hz 
in several of our calculations. We flag such cases when 
appropriate. 

To understand more clearly how LISA extracts the lu- 
minosity distance from measurements of a binary black 
hole inspiral, it is useful to rewrite somewhat schemati- 
cally the measured form of the inspiral as follows: 

M 5/3 f(t) 2 / 3 

h?u%t) = z ^,„( "angles", t) cos[$(t)+^, n ( "angles", t)\ 

We have subsumed the angle-dependent factors (L ■ n) 
and F^i* into the schematic functions !F\n{ "angles", t); 
we leave the dependence upon t in these functions as a 
reminder that the constellation's motion modulates the 
waveform. We have likewise written the phase modula- 
tions imposed by the detector's response and motion in 
the schematic form ipi,n ( "angles", t). 

From the form of equation JSJ we see that the lumi- 
nosity distance is very strongly correlated with the red- 
shifted chirp mass, M z , and the various angles which set 
the instantaneous waveform amplitude. As already men- 
tioned, M z is typically determined with extremely high 
precision because it fixes the phase evo l ution: typically , 
SM Z /M Z < 0.01%. See lHughesI (12001 . IVeTx-hiol ll200l 
for examples specific to LISA. 

The modulations induced by LISA's orbital motion 
make it possible to measure sky position for events which 
last for at least a fair fraction of LISA's orbit. We es- 
timate the accuracy with which position (among other 
parameters) is determined using a maximum likelihood 
parameter estimation formalism l)Finnl ll992): from the 
detector's response to a given gravitational wave, we con- 
struct the variance-covariance matrix E ah . Diagonal ele- 
ments of this matrix represent the rms error ({Sx a ) 2 ) in 
a source parameter % a ; off-diagonal components describe 
the degree to which errors in parameters x° & n d X b are 
correlated. See iHugriesI l)2002t) for discussion specifically 
tailored to this application. 

2.3. Measurement accuracy distributions 

To assess distance and position accuracy, we have es- 
timated the accuracy with which these parameters are 
measured for a wide range of binary masses. For each 
set of masses, we randomly distribute the sky position 
and orientation of 10,000 such binaries. We then calcu- 
late the fractional accuracy with which distance is deter- 
mined for each binary, SDl/Dl, as well as the angular 
sky position error 69. Figure [21 shows the distribution we 
find in these quantities for binaries with mi = 10 5 Mq, 
m 2 = 6 x 10 5 M at z = 1. We find that the typical po- 
sition determination is relatively poor — these binaries 
are fixed to an error box that, at best, is ~ 5 arcminutes 
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m, = 10 B M_, m„ = 6xlO B M. at z = 1 



m, = 10 5 M_, m, = 6xl0 5 M n at z = 3 
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Fig. 2. — Pointing and distance error distributions for mea- 
surements at z = 1 of a binary of masses mi = 10 s M©, mi = 
6 X 10 5 Mq . These distributions were made by Monte-Carlo simu- 
lations of 10,000 LISA BBH measurements, randomly distributing 
the bi naries' positions, orientations, and merger times; see Huehes 
12002) for details. The top distribution shows that the most likely 
position error boxes have sides 58 < 10 arcminutes, spreading out 
to 58 > 3°. The distance distribution peaks at 5Dl/Dl < 1%, 
with most of the distribution confined to SD^/D^ < 5%. 




Fig. 3. — Pointing and distance error distributions for mea- 
surements at z = 3 of a binary of masses mi = 10 s Mq , 7712 = 
6 X 10 s Mq. The distribution for position error is so broad that 
we cannot really identify a "most likely" position error; however, 
most of the distribution lies at 88 < 10° . The distance distribution 
peaks at SD^/D^ < 10%, with most of the distribution confined 
to 5D L /D L <30%. 



on a side. In most cases, the resolution is substantially 
worse. The distance determination, by contrast, is quite 
good: half of these events have their distance determined 
with precision SDl/Dl < 1%. 

Table 1 summarizes the parameter determination dis- 
tributions we find for a wide range of masses. For 
SD L /D L and 59, we give the 5%, 25%, 50%, and 90% 
likelihood values from the distributions predicted by our 
Monte-Carlo calculation. For example, for binaries with 
mi = 3 x 10 4 M Q , m 2 = 10 5 M at z = 1, 25% of all 
events have SDl/Dl < 0.006 and localize the source 
to 59 < 14.6 arcminutes; 90% of all events with these 
masses and redshifts have 5Dl/Dl < 0.029 and localize 
the source to 59 < 120 arcminutes. 

Notice that, in this table, the best pointing and dis- 
tance determination occurs for binaries that have a total 
(redshifted) mass (1 + z)(mi + 7712) ~ several x 10 5 Mq. 
Two competing effects drive this behavior. First, for 
small binaries the amplitude of the GWs is smaller; the 
degradation of their parameter determination is due to 
reduced signal-to-noise ratio. Larger binaries are more 
interesting. When such binaries enter LISA's sensitive 
band, they are closer to their final merger — much 
less inspiral remains once LISA begins measuring their 
waves. They therefore do not exhibit as many cycles 
of detector-motion-induced modulation, and so their po- 
sition angles are not as well determined. In particu- 
lar, we find that distance and position determination 
rapidly degrades as binaries are made more massive than 
(1 + z)(mi +m 2 ) > a few x 1O 6 M . 

The poor parameter determination tendency of large 
binaries can be repaired somewhat by improving LISA's 
low frequency sensitivity. If the antenna has good sen- 
sitivity at lower frequencies, the span of data containing 



good information about the inspiral can be lengthened. 
Table 2 shows how well we measure distance and position 
when fi ow is reduced from 10~ 4 Hz to 3 x 10 -5 Hz. We 
now find that the distance is determined very precisely 
for binaries with total mass (several) x 10 6 Mq. Sky posi- 
tion error is no worse than that achieved at lower masses 
- 59 < 10 arcminutes in the best cases, and is more 
typically a factor of a few larger than this. 

The same general story holds as we move to larger 
rcdshift. Figure duplicates the content of Figure[3 but 
with binaries at redshift z = 3. Likewise, Tables 3 and 
4 duplicate the content of Tables 1 and 2, respectively, 
with all binaries placed at z = 3. The overall parame- 
ter determinations are worsened, as we would expect — 
these sources are much farther away, and so have greatly 
reduced signal-to-noise. In addition, the larger cosmolog- 
ical redshift shifts the signal to lower frequencies, where 
much of it is lost in low-frequency noise. To quantify the 
impact of this effect, in Table 4 we present results show- 
ing what happens when we lower /i ow from 10 -4 Hz to 
3x 10~ 5 Hz. All cases with mi+m 2 > 6x 10 5 M Q are sub- 
stantially improved by this fix. Good low frequency per- 
formance will be important for measuring high redshift 
binaries. The best pointing accuracy we find is 59 < 40 
arcminutes; 59 ~ 1° or larger is more typical. The lu- 
minosity distance can still be determined quite well - 
we find errors of a few percent in the best cases, and 
5Dl/Dl < 15% is quite common. 

It is worth emphasizing at this point that the re- 
sults we present here are most likely somewhat conser- 
yative. By taking into account other GW harmonics 
ijMoore fc Hellingsl2002|) and properly accou nting for th e 
high-frequency structure of LISA's response (|Setol l2002'l . 
the pointing accuracy, and thus distance accuracy, can 
be improved by a factor of a few. Properly accounting 
for modulations induced by spin-orbit and spin-spin cou- 
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Fig. 4. — Distance errors for BBH measurements at z = 1 with 
mi = 10 5 Mq, ni2 = 6x 10 5 Mq, assuming that an electromagnetic 
counterpart allows precise sky position determination. The peak 
error is at SD^/D^ ~ 0.1%, and is almost entirely confined to 
5D L /D L <0.5%. 
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Fig. 5. — Distance errors for BBH measurements at z = 3 with 
mi = 10 5 Mq, m2 = 6x 10 5 Mq, assuming that an electromagnetic 
counterpart allows precise sky position determination. The peak 
error is at SD^/D^ ~ 0.5%, and is almost entirely confined to 
SD L /D L <2%. 



pling can also improve pointing accuracy and t hus dis- 
tance determination, in some cases significantly ijVecchiol 
l2f)0l . 

Using a determination of Dl , we can infer the redshift 
by using knowledge of the universe's geometry (the Hub- 
ble constant, mean density of mat ter O m , and density 
of dark energy Clx) l)Hughesll2002j) . This makes possi- 
ble interesting analyses (e.g., we can map the distribu- 
tion of black hole masses as a function of redshift) , but 
presupposes rather than measures the distance-redshift 
relation. 

3. THE IMPACT OF A COUNTERPART 

Parameter estimation improves dramatically when an 
EM counterpart to a BBH GW event can be identified. 
The counterpart will almost certainly be pinpointed with 
far greater accuracy than is possible with GWs. Corre- 
lations between position and distance are then broken, 
greatly reducing the distance error. An example of this 
improvement is shown for z = 1 in Figure El The dis- 
tribution of distance errors peaks near SDl/Dl ~ 0.1%, 
and is largely confined to SDl/Dl < 1%. Similar results 
are seen for z = 3 (Fig. [5J , albeit with precision degraded 
by a factor of a few due to lower signal-to-noise. Com- 
paring to the lower panel of Figures El and El we see that 
associating the event with a counterpart improves dis- 
tance accuracy by roughly an order of magnitude. This 
rough level of improvement holds over a wide band of 
mass and redshift. 

A correlated GW/EM measurement will be particu- 
larly important if the counterpart provides a redshift as 
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Fig. 6. — Likelihood contours for measurement of the matter 
density H m and dark energy equation of state parameter w (with 
the pressure and density of the dark energy related by p = wp). 
We assume that the universe is flat, and that the underlying model 
has f2 m = 0.3 and w = — 1. The two GW sources are at z = 1 
and z = 3, while the SNAP SNe are evenly distributed within 
0.7 < z < 1.7. 



well as an improved Dl ■ Such a measurement would con- 
stitute a powerful standard candle, probing the distance- 
redshift relation in a manner complementary to other 
candles, such as Type la supernovae. Assuming a flat 
universe and a Hubble constant h$ = 0.65, we ask how 
well the matter density £l m and dark energy equation-of- 
state ratio w can be measured. Figure El shows estimated 
likelihood contours (1-cr) in the (Q m ,w) plane. The dot- 
ted line shows the contour expected for measurements 
of 3,000 SNe evenly distributed within 0.7 < z < 1.7 
(reasonable choices for SNAP); the solid contour is for 
two GW events, one each at z = 1 and z = 3. (We 
will discuss the dashed line further below.) Redshift and 
distance are measured with such accuracy that the con- 
tours are extremely tight even for only a small number 
of sources. 

We emphasize the current poor understanding of EM 
counterparts to BBH GW events, although the possi- 
bility of such counterparts has been discussed for quite 
some time (e.g., Begelman, B landford, & Rees 1980). 
iMilosavlievic & Phinnevl l)2004() have recently examined 
the evolution of gas in the environment of a merging bi- 
nary, and show that there is likely to be a delayed electro- 
magnetic afterglow. They find that the merging binary 
carves a hollow region in the volume of circumbinary 
gas. The binary separation shrinks faster than the in- 
ner edge of the hollowed region; thus, as the coalescence 
proceeds, there should be no substantial accretion of ma- 
terial onto the system. The gas falls onto the merged 
remnant several years after the merger, leading to an af- 
terglow that should be measurable by next generation 
x-ray telescopes. 

Other models suggest that there may be an electromag- 
netic precursor to the merger, rat her than a delayed glow. 
One e xample is discussed by l Armitage fc Nataraianl 
(2002). They argue that gas is driven onto the larger 
member of the binary by the secondary's inspiral, lead- 
ing to super-Eddington accretion. In this model, much of 
the inner disk may be expelled from the system in a high 
velocity (~ 10 4 km/sec) outflow. Such strong outflows 
could flag a recent or impending merger. A similar family 
of models (|Sillanpaa et alJl988tlLehto fc Valtonenl 199611 
explains periodic variations in the BL Lac object OJ 287 
by a tight, eccentric binary system with mass ratio of 
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Fig. 7. — The differential probability of magnification by grav- 
itational len sing, p(ai), for sources at z = 1.5 in a concordance 
universe; see Wang, Holz, & Munshi (2002) for details. 



about 1:100. Flaring outbursts from this quasar are ex- 
plained as arising from the secondary periodic crossing 
of the primary's accretion disk. Given the great payoff 
that would follow from associating a counterpart to a 
GW event, we strongly advocate continuing to develop 
and refine models of BBH mergers. 

It is worth noting that, for a small fraction of bina- 
ries (assuming a sufficiently high event rate), LISA will 
provide an error box of < 5 arcmin, and an estimate of 
the time of merger about a day in advance. Regardless 
of the state of theoretical predictions, we imagine that 
in such cases there will be great interest in searching the 
GW source error box for any observational counterparts 
to the merger. Indeed, as we briefly discuss in Sj5j the 
number of relevant galaxies in the LISA error box may 
be fairly small, and so associating an EM counterpart to 
the GW event may be tractable. 

4. GRAVITATIONAL LENSING 

Having discussed the impressive quality of GW stan- 
dard sirens, we turn now to an important caveat: the 
impact of gravitational lensing on the distance measure- 
ment. GWs are lensed exactly as EM radiation is lensed. 
Since we expect BBH events to come from rather large 
rcdshift [z > 1), weak lensing in the GW datasets should 
be common l|Markovi6!ll993t iWang. Stebbins. fc Turned 
1996) (in addition to the occasional strongly-lensed 
source) . 

A lens with magnification /i will distort the inferred 
luminosity distance to the source: if the true distance 
is Dli we measure D^/y/Ji, incurring a "systematic" er- 
ror ADl/Dl = 1 — 1/v 7 ^- We estimate the error such 
lensing is likely to introduce by convolving this quan- 
tity with the expected magnification distribution, p(fi) 
l|Holz fc Waldl Il998t IWang. Holz. fc Munshil jflOg) ; an 
example of this distribution is shown in Figure [7| Us- 
ing parameters appropriate to a ACDM model of the 
universe, we find a mean error at z = 2 of (ADl/Dl) — 
0.005, with a variance y/ ((AD L / D L ) 2 ) ~ 0.09 (approx- 
imatin g the lensing by a Gaussian; see Ho lz fc Linderl 
( 2005)). The dashed line in Figure El shows the contour 
we expect from the 2 GW sources when lensing errors 
are included. The parameter accuracies are significantly 
degraded. 

Of course, this magnification bias affects all standard 
candles, not just GWs. The rate of Type la SNe, how- 
ever, is high enough to sufficiently sample the entire lens- 
ing distribution, and thus average away the bias. Mis- 
sions such as SNAP are designed to observe thousands 
of SNe at high redshift, in large part to overcome gravi- 



tational lensing. Indeed, this may allow one to measure 
the lensing signal w ell enough to infer characteristics of 
the l ensing matter (jMetcalf fc Silklll999t iSeliak fc HoTzl 
|1999|) . This is unlikely to be the case with BBH GWs: 
the r ate of mergers will likely be muc h lower than that of 
SNe l|Rlchstonell998HHaehne1Hl 9981 , so we cannot count 
on enormous numbers of events. We also emphasize that 
we do not expect to be able to corre ct for gravitationa l 
lensing effects on a case- by-case basis (Dalai et al. 2003). 
Lensing, therefore, will introduce an insurmountable er- 
ror of ~ 5-10% for each individual high redshift event, 
significantly greater than the intrinsic distance error. 

5. IDENTIFYING THE COUNTERPART 

In order to provide data on the distance-redshift curve, 
a GW event must be associated with an "electromag- 
netic" counterpart — GWs provide an accurate measure 
of luminosity distance, but give no direct information 
about redshift. This is the weakest link in our analy- 
sis: we do not know whether such counterparts exist. 
However, a simple counting argument suggests that the 
number of relevant galaxies in the LISA error cube may 
be fairly small. We approximate the redshift distribution 
of source galaxies by 



dN 



cx iTexp [~(R/R*f] , 



(6) 



dRdtt 

where R is the co moving distance; we ta ke a = 1, /3 = 4, 
and R* = c/H l|Kaiserlll992l lHullT999|) . We normalize 
this to a projected number density of 



dN 

In 



dR 



dN 
dRdVL 



300 galaxies/ arcmin 



(7) 



appro ximating the Hubble Deep Field fWilliams. et alJ 
1996). 

As we have discussed extensively, a GW measurement 
of a binary black hole merger determines the position on 
the sky to within some error 59, and determines the lumi- 
nosity distance to within some error SD^. By assuming 
a cosmological model we can convert the measured lumi- 
nosity distance, and its error, to any other desired cosmic 
distance measure. Denoting by 5R(5Dl; <5cosmology) 
the error in comoving distance due to both the GW 
measurement error and the uncertainty in cosmological 
parameters, the number of galaxies which lie in the 3- 
dimcnsional GW error cube is 



N 



dN 



error cube 



, 56 2 SR(6D L ; Cosmology) . (8) 

dRdil 

Figure IS] shows four realizations of N eTTOI cu b c as a func- 
tion of GW event redshift, z. We have scaled to a near 
best-case pointing accuracy of 59 = 1 arcmin; we em- 
phasize that this is an optimistic, though not implausi- 
ble, pointing error. From the Tables it is apparent than 
about 5% of binaries with masses and redshifts in the 
LISA sweetspot have positional errors 59 < 5 arcmin. As 
discussed towards the end of H2.3I these results may be 
conservative — accounting for spin-induced precessional 
effects may allow certain degeneraci es to be broken and 
improve LISA's pointing accuracy lIVecchiolEiol . [A 
fiducial pointing error of 1 arcminute also makes it very 
simple, by Eq. ©, to scale to larger values.] The number 
of galaxies decreases significantly as uncertainties in cos- 
mological parameters are reduced (as is to be expected by 



Using gravitational-wave standard sirens 



7 



current + lensinc 
current 
future + lensing 
future 




Fig. 8. — Number of candidate host galaxies per square ar- 
cmin of LISA error cube for a supermassive binary black hole co- 
alescence event, as a function of redshift of the event. Dotted 
(gree n) line utilizes curr ent uncertainties in cosmological parame- 
ters ISperge l ct al. 2003). Dash-dotted (blue) line represents pos- 
sible future improvements in these parameters (1% in f! m , 
ho). Dashed (red) line includes the degradation in the depth of the 
error cube due to gravitational lensing, for current cosmological 
uncertainties. Solid (black) line represents future cosmological un- 
certainties, with the inclusion of gravitational lensing degradation. 
The 1 arcmin LISA sky position error considered here is optimistic, 
though not implausible (see text). This choice allows for a straight- 
forward scaling of the curves to larger positional errors. 



the time that LISA is operating). In Figure[S]we include 
the curve for the current state of errors on cosmological 
parameters, as well as for expected future, percent-level, 
measurements. As was discussed in the previous section, 
gravitational lensing adds a further, insurmountable er- 
ror to the GW measurement of SD^. We approximate 
the lensing effects by a Gaussian i n magnification, wit h 
variance given by censing = 0.088z (|Holz fc Linderl2 005). 
Although this expression is strictly appropriate only for 
high source statistics (as otherwise the lensing distribu- 
tion is non-Gaussian), it is a sufficiently good approxima- 
tion for present purposes. With the inclusion of lensing 
errors in addition to cosmological parameter uncertain- 
ties, we find < 10-20 potential counterpart galaxies in 1 
arcmin of LISA error cube. 

The number of candidate host objects for a galactic 
binary black hole merger is thus likely to be tractable. 
What remains is to find a way to identify which of 
the dozen or so candidate objects is in fact the host 
of the merger. Useful models exploring the signatures 
that may constitute "precursor s " of the merger ex- 
ist llArmitaee fc Nataraianll2002l iSillannaa et alJ 119881 
ILehto fc Valtonenl ll996) ; we are hopeful that such mod- 
els can be extended to the ranges of mass and mass ratio 
of binaries that are likely to be observationally interest- 
ing to LISA. 

The remnant of the merger is very likely to have an ir- 
regular morpholog y. In addition, the host ga l axy m ay be 
in an active phase. Milosavl ievic fc Phinnevl ((2004) have 
recently develop ed a model fo r the la te x-ray afterglow of 
a BBH merger. IKocsis et alJ l|2005j) consider a scenario 
in which a merger remnant is associated with a quasar, 
and argue that in this case the paucity of quasars will 
make the identification of a counterpart significantly eas- 
ier. In both of these models the merger remnant "lights 
up" , and is thus relatively easy to identify in the posi- 
tional error box. Even if the remnant does not light up, 



other information will help winnow the list of candidate 
host galaxies. For example, the GWs will measure the 
masses of the black holes with good accuracy; by using 
properties such as the Mbh~c relation we can estimate 
specific properties (e.g., kinematics, luminosity) that the 
host galaxy would be expected to have. 

Observations of the LISA error boxes will no doubt be 
undertaken, regardless of the state of theoretical predic- 
tions. It is only by such direct observations that we will 
determine whether or not EM counterparts to BBH GW 
sources can be identified. 

6. CONCLUSIONS 

Because of their potential as an independent set of 
standard candles, BBH GW standard sirens can make an 
important contribution to programs to map the distance- 
redshift relation over a large span of redshift. Although 
the intrinsic precision of these candles is phenomenal, 
this precision will be limited in practice because of grav- 
itational lensing. With lensing taken into account, the 
accuracy of the BBH GW candle is comparable to (or 
perhaps slightly better than) a Type la SN. It is sober- 
ing to note that we are already approaching the point at 
which lensing, rather than intrinsic dispersion, limits our 
ability to use standard candles. 

We emphasize that the systematics of BBH events are 
entirely different from those of Type la SNe. As such, 
the greatest impact of BBH standard sirens may be to 
verify, and thereby increase our confidence in, other stan- 
dard candles. The utility of these GW standard sirens 
depends on the identification of an electromagnetic coun- 
terpart, through which a redshift to the source can be 
determined. It is not unlikely that, at least for some of 
the best-observed systems, a counterpart will be found. 
If this is the case, the BBH GW source would become an 
exceptionally precise standard siren. 
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TABLE 1 

Measurement precision at z = 1 with /i ow = 10 -4 Hz 



mi (M ) 


m 2 (M ) 


SD L /D L 


89 








(5%, 25%, 50%, 90%) 


(5%, 25%, 50%, 90%) 


10 4 


10 4 


(0.005, 0.010, 0.014, 0.042) 


(14.3, 28.3, 48.5, 117) 


arcmin 


10 4 


3 x 10 4 


(0.003, 0.008, 0.013, 0.037) 


(11.0, 22.1, 41.6, 111) 


arcmin 


10 4 


6 x 10 4 


(0.003, 0.007, 0.013, 0.036) 


(9.05, 19.1, 40.4, 109) 


arcmin 


10 4 


10 5 


(0.005, 0.007, 0.013, 0.035) 


(7.85, 18.3, 39.5, 110) 


arcmin 


3 x 10 4 


10 5 


(0.002, 0.006, 0.013, 0.029) 


(5.26, 14.6, 37.4, 120) 


arcmin 


6 x 10 4 


10 5 


(0.002, 0.006, 0.012, 0.037) 


(4.07, 12.9, 35.0, 120) 


arcmin 


10 5 


10 5 


(0.002, 0.005, 0.012, 0.034) 


(3.36, 11.7, 33.4, 117) 


arcmin 


10 5 


3 x 10 B 


(0.001, 0.005, 0.012, 0.035) 


(2.71, 10.6, 33.1, 116) 


arcmin 


10 5 


6 x 10 5 


(0.001, 0.006, 0.014, 0.044) 


(2.82, 12.0, 38.8, 120) 


arcmin 


10 5 


10 6 


(0.002, 0.009, 0.017, 0.053) 


(3.89, 18.5, 50.9, 126) 


arcmin 


3 x 10 5 


10 6 


(0.002, 0.013, 0.026, 0.087) 


(4.65, 29.7, 71.0, 172) 


arcmin 


6 x 10 5 


10 6 


(0.003, 0.019, 0.035, 0.122) 


(5.60, 39.2, 93.6, 220) 


arcmin 


10 6 


10 6 


(0.004, 0.024, 0.043, 0.149) 


(6.36, 52.2, 118, 271) 


arcmin 



TABLE 2 

Measurement precision at z = 1 with / 1ow = 3 x 10 -5 Hz 



mi (M Q ) 


m 2 (M Q ) 


SD L /D L 


56 








(5%, 25%, 50%, 90%) 


(5%, 25%, 50%, 90%) 


10 5 


10 5 


(0.002, 0.005, 0.012, 0.036) 


(3.35, 11.4, 33.2, 117) 


arcmin 


10 5 


3 x 10 5 


(0.001, 0.005, 0.011, 0.034) 


(2.68, 10.4, 31.0, 107) 


arcmin 


10 5 


6 x 10 5 


(0.001, 0.005, 0.010, 0.032) 


(2.71, 9.85, 29.7, 103) 


arcmin 


10 5 


10 6 


(0.001, 0.005, 0.010, 0.031) 


(3.27, 10.5, 29.1, 101) 


arcmin 


3 x 10 5 


10 6 


(0.001, 0.004, 0.009, 0.028) 


(2.45, 8.77, 25.8, 89.9) 


arcmin 


6 x 10 5 


10 6 


(0.001, 0.004, 0.009, 0.026) 


(2.25, 8.5, 24.8, 84.3) 


arcmin 


10 6 


10 6 


(0.001, 0.004, 0.009, 0.026) 


(2.39, 9.05, 25.3, 83.5) 


arcmin 


10 6 


3 x 10 6 


(0.002, 0.006, 0.011, 0.034) 


(4.56, 12.6, 34.6, 93.8) 


arcmin 


10 6 


6 x 10 6 


(0.003, 0.008, 0.016, 0.051) 


(7.16, 20.6, 45.9, 110) 


arcmin 


10 6 


10 7 


(0.003, 0.011, 0.021, 0.071) 


(9.20, 27.0, 57.1, 136) 


arcmin 



TABLE 3 

Measurement precision at z = 3 with / 1ow = 10 -4 Hz 



mi (Mq) 


m 2 (M e ) 


SD L /D L 


89 








(5%, 25%, 50%, 90%) 


(5%, 25%, 50%, 90%) 


10 4 


10 4 


(0.013, 0.029, 0.051, 0.143) 


(37.7, 78.2, 158, 422) 


arcmin 


10 4 


3 x 10 4 


(0.010, 0.026, 0.050, 0.135) 


(27.2, 66.8, 150, 428) 


arcmin 


10 4 


6 x 10 4 


(0.008, 0.024, 0.050, 0.145) 


(22.0, 57.7, 140, 442) 


arcmin 


10 4 


10 5 


(0.008, 0.024, 0.050, 0.142) 


(19.3, 55.4, 143, 465) 


arcmin 


3 x 10 4 


10 5 


(0.006, 0.019, 0.044, 0.131) 


(12.6, 44.5, 125, 444) 


arcmin 


6 x 10 4 


10 5 


(0.006, 0.021, 0.044, 0.128) 


(10.1, 42.6, 132, 429) 


arcmin 


10 5 


10 s 


(0.005, 0.024, 0.049, 0.141) 


(9.55, 44.8, 144, 430) 


arcmin 


10 5 


3 x 10 s 


(0.007, 0.034, 0.069, 0.213) 


(11.3, 69.9, 193, 485) 


arcmin 


10 5 


6 x 10 s 


(0.008, 0.044, 0.087, 0.287) 


(17.5, 96.1, 240, 593) 


arcmin 


10 5 


10 6 


(0.009, 0.058, 0.111, 0.378) 


(23.4, 127, 304, 734) 


arcmin 
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TABLE 4 

Measurement precision at z = 3 with /i ow = 3 x 10~ 5 Hz 



mi(M ) m 2 (M Q ) 5D L /D L 69 

(5%, 25%, 50%, 90%) (5%, 25%, 50%, 90%) 





10 4 




10 4 


(0.013, 


0.029, 


0.050, 


0.136) 


(37.4, 76.0, 156, 423) 


arcmin 




10 4 


3 


x 10 4 


(0.010^ 


0.026^ 


0.050^ 


0.137) 


(26.5^ 62. 7, U9, 437) 


arcmin 




10 4 


6 


x 10 4 


(0.009, 


0.024, 


0.049, 


0.140) 


(21.9, 57.8, 141, 441) 


arcmin 




10 4 




10 5 


(0.007, 


0.023, 


0.049, 


0.139) 


(19.0, 55.3, 143, 465) 


arcmin 


3 


x 10 4 




10 5 


(0.006, 


0.019, 


0.043, 


0.135) 


(12.4, 43.6, 124, 429) 


arcmin 


6 


x 10 4 




10 5 


(0.005, 


0.018, 


0.040, 


0.126) 


(10.4, 38.8, 114, 404) 


arcmin 




10 5 




10 5 


(0.004, 


0.016, 


0.038, 


0.120) 


(8.85, 34.7, 106, 384) 


arcmin 




10 5 


3 


x 10 5 


(0.004, 


0.016, 


0.035, 


0.108) 


(8.26, 32.4, 98.7, 360) 


arcmin 




10 5 


6 


x 10 5 


(0.005, 


0.016, 


0.035, 


0.105) 


(10.9, 36.0, 100, 346) 


arcmin 




10 5 




10 6 


(0.006, 


0.017, 


0.035, 


0.106) 


(15.2, 41.6, 105, 349) 


arcmin 


3 


x 10 5 




10 6 


(0.005, 


0.017, 


0.034, 


0.100) 


(12.2, 37.3, 101, 327) 


arcmin 


6 


x 10 5 




10 6 


(0.006, 


0.019, 


0.039, 


0.115) 


(13.1, 42.1, 120, 336) 


arcmin 




10 6 




10 6 


(0.007, 


0.025, 


0.049, 


0.146) 


(16.9, 55.5, 145, 359) 


arcmin 




10 6 


3 


x 10 6 


(0.009, 


0.040, 


0.077, 


0.240) 


(24.8, 89.1, 208, 507) 


arcmin 




10 6 


6 


x 10 6 


(0.014, 


0.063, 


0.116, 


0.392) 


(35.4, 142, 322, 757) 


arcmin 



